Observation of Josephson coupling through an interlayer of antiferromagnetically 

ordered chromium 
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The supercurrent transport in metallic Josephson tunnel junctions with an additional interlayer 
made up by chromium, being an itinerant antiferromagnet, was studied. Uniform Josephson coupling 
was observed as a function of the magnetic field. The supercurrent shows a weak dependence on the 
interlayer thickness for thin chromium layers and decays exponentially for thicker films. The diffusion 
constant and the coherence length in the antiferromagnet were estimated. The antiferromagnetic 
state of the barrier was indirectly verified using reference samples. Our results are compared to 
macroscopic and microscopic models. 

PACS numbers: 



I. INTRODUCTION 

The field of superconducting spintronics comprises in- 
teresting physical phenomena with potential apphcations 
for digital and quantum logics. The Cooper pairs leaking 
from a conventional superconductor (S) into a ferromag- 
net (F) display phase oscillation of their order parameter 
and a rapid decay of the amplitude over a few nanome- 
ters inside the F-laycr [1]. These phase oscillations are 
for example used to construct tt coupled S-F-S Joseph- 
son junctions (JJs), for which the Josephson phase in 
the ground state is tt instead of 0, or O-tt coupled JJs, 
where the and tt coupling is locally set by a stepped 
F-layer [2]. If the supercurrent leaks into an itinerant 
antiferromagnet (AF), the spin-dependent quasiparticle 
and Andrecv reflections were shown to create low-energy 
bound states [ '] leading for example to atomic-scale 
to TT transitions of the coupling in S-AF-S type JJs [4]. 
S/AF multilayers are model systems for antiferromag- 
netic superconductors because pairing and pair breaking 
can be locally separated. 

The to TT phase oscillation in S-F-S JJs was veri- 
fied in a number of publications [■"), G, 7], but up to now 
only a few experiments have been performed on S-AF-S 
JJs. Fast oscillations of the critical current Ic versus the 
magnetic field H were observed over an oxide AF inter- 
layer and explained by canting of its magnetic moments 
(similar to the giant magnetoresistance effect) [8]. Trans- 
port studies in metallic S-AF-S films employing a disor- 
dered FeMn alloy as AF [9] showed deviations from the 
conventional Ic{H) pattern, too, indicating (i) a nonuni- 
form current distribution or (ii) a local change in the 
AF magnetization. Transport measurements on S-AF-S 
junctions based on mono-atomic chromium or its alloys 
were proposed in Refs. 4, 10, but not realized yet. 
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FIG. 1: Schematic representation of a S-AF-S type Josephson 
junction with longitudinal spin density waves on the two anti- 
ferromagnetically coupled sublattices of the Cr interlayer (low 
temperature state), for which both the propagation vector Q 
and the spins point out-of-plane. 



In this paper we study the supercurrent transport 
through an antiferromagnetically ordered interlayer. The 
3d metal chromium as one of the three elemental anti- 
ferromagnets (Cr, a-Mn, 7-Fe) apart from the rare earth 
and actinoids elements, was chosen due to its simple crys- 
talline structure and low constrains on its atomic order 
for the onset of antiferromagnetic order. 

The magnetic field dependence of the critical current 
indicates uniform Josephson coupling over the junction 
area for the studied magnetic field range up to 0.5 mT. 
For thin AF thicknesses the supercurrent was found to 
scatter strongly but to be only weakly dependent on the 
chromium thickness, whereas for thicker interlayer thick- 
nesses the supercurrent decays exponentially. The super- 
current coherence length, diffusion constant, and mean 
free path were determined. The magnetic properties of 
the Cr interlaycrs were verified by SQUID magnetometry 
using reference samples. 
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II. THEORY 

Itinerant antiferromagnets like Cr with spin density 
waves (SDW) propagation along special crystal direc- 
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tions [1 y] and nesting properties are described by BCS- 
like triplets of electron-hole pairing. Two pieces of the 
Fermi surface are connected by the SDW wavevcctor Q, 
and the electron and hole surfaces can be superposed by 
translation about Q (so called nesting condition). The 
orientation of SDW depends strongly on the tempera- 
ture and sample properties like the magnetic and atomic 
interface structure [1 ) ]. In the low temperature limit and 
for a clean interface to a FM the SDW in Cr is longitu- 
dinal, i.e. its spins and wavevector are oriented parallel 
and point out-of-plane, see Fig. 1 [12]. Note that for the 
transport of Cooper pairs the orientation of the quanti- 
zation axis is not of importance, and may be orientated 
differently in the antiferromagnetic interlayers of our JJs. 
The incommensurable spin density oscillation length of 
Cr at low temperature (< 80 K) is about 42 monolayers 
or 6 nm [11]. Internal elastic strain or small grain sizes, 
as inevitably present in polycrystalline films, may cause 
the SDW to become commensurate. Bulk Cr has an AF 
exchange energy Ee^ = 120 meV, a Neel temperature 
T/v = 311 K, and an atomic magnetic moment of 0.5 /ie 
at 4.2 K. The AF exchange energy E^j; is much larger 
than the superconducting gap energy, and the mutually 
compensating magnetic moments prevent spin accumu- 
lation in the superconductor when biasing S-AF-S JJs. 

The macroscopic model, based on the quasiclassical 
theory of an antiferromagnetic interlayer with nesting 
condition, like Cr, in the dirty limit £ <C £,af (mean free 
path £, coherence length £,af) and close to the critical 
temperature of the superconductor T^, is described by 
the linearized Usadel equation [l-;!] with the anomalous 
Green's function Faf inside the AF-layer [10]: 

_hDAF p ... 

with an exponentially decaying solution of the form 

Faf ^ exp {-cLafHaf) (2) 

and coherence length ^af = \/^J^af/^^E^ with diffu- 
sion constant Daf = vpi/^, mean free path i and AF- 
layer thickness (Iaf- This ansatz is similar to the con- 
ventional form for ferromagnets [l] and differs by a fac- 
tor 2 from the solution used in Refs. 9, 10. By replac- 
ing the AF-layer with a F-layer the magnetic exchange 
field Eex in Eq. (1) becomes imaginary and its solu- 
tion exp ^—dp/ i^p + contains decay and oscilla- 
tions lengths = ~ \/ hD p / E^x {dp- F-layer thick- 
ness) . The decay length is by a factor smaller than 
S,AF in a AF-layer [10]. The solution for Faf. Eq. (2), 
cannot provide a change of the Josephson ground state 
from to TT phase. 

The microscopic model, taking the atomic magnetic 
order of the AF-layer into account, describes the ground 
state phase (0 or tt) by the spin-up and spin-down re- 
flection amplitudes. For an odd number of AF-layers tt 
coupling can be obtained, whereas for an even number of 
layers the JJs are always in the ground state [4]. 



In Josephson junctions with an additional tunnel bar- 
rier (I), i.e. SI-AF-S JJs as studied in this paper, the 
low transparency interface on one side of the AF inter- 
layer modifies the density of states profile in the super- 
conductor and the Andreev reflections in the interlayer. 
For example, in the case of JJs with a ferromagnetic in- 
terlayer the /c(r, c^f) dependencies for S-F-S differ from 
SI-F-S JJs [14]. However, to our knowledge, neither (i) 
S-AF-S junctions with one low transparency interface, 
i.e. SI-AF-S JJs nor (ii) itinerant AFs with incommen- 
surable SDWs are explicatively included in the currently 
available theoretical models. 



III. EXPERIMENT AND DISCUSSION 

The sputtered SI-AF-S multilayers consist of Nb elec- 
trodes and an AlOx tunnel barrier (I). The AF layer was 
wedge-shaped, thus all JJs of one set were fabricated and 
patterned simultaneously [15]. A 2nm Cu buffer-layer on 
top the I-layer improved the growth of the antiferromag- 
netic Cr due to the very low solubility of Cr in Cu [hi]. 
Thus, the stack was SIN-AF-S-type hke, with the thin N- 
layer not affecting the junction parameters as determined 
from SIS and SINS reference samples [15]. The AlOa: 
tunnel barrier was formed dynamically at 1 • 10"'^ mbar 
(sample sets 1 and 3) or 3 • 10~^ mbar (set 2) residual 
oxygen pressure. Similar SI-F-S-type JJs showed a uni- 
form increase of the average interlayer thickness and low 
interface interdiffusion despite the polycrystalline growth 
[6, 15, 17, 18]. The lithographically patterned JJs had 
areas of 10 x 5 and 50 x 10 /im^, and an effective length 
ranging from the intermediate to the short JJ limit, i.e. 
L/A,7 = [4 . . . 0.1] with the Josephson penetration length 
A J ~ I/a/Tc- Inserting a tunnel barrier in the S-AF-S 
stack increases the normal state and subgap resistances. 
Thus, DC measurements of these samples are more fea- 
sible and the Josephson dynamics can be observed, as 
seen by the underdamped /^-characteristics in the inset 
of Fig. 2. 

The SI-AF-S JJs were zero-field cooled down to 4.2 K 
using /i-metal shields to suppress the external stray fields. 
The magnetic field H was applied in-plane and the bias 
current was computer-controlled swept while measuring 
the voltage drop across the junction. Room-temperature 
voltage amplifiers were used. Both current and voltage 
values were automatically averaged over several hundred 
data points. The resolution limit is = 1 /iA. 



A. Current-voltage and magnetic field dependence 

The STAF-S JJs had hystcretic current-voltage char- 
acteristics up to dAF ^ 15 nm with reduced subgap fea- 
tures compared to normal SIS-type JJs, see inset of Fig. 
2. The position of zero-field and Fiske steps were as ex- 
pected and were not further investigated in this work. 
The magnetic field dependence of the Fraunhofcrn pat- 
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FIG. 2: (Color online). Normalized Ic(H) curves and IV- 
characteristics (inset) for various (Iaf from sets 1 and 2. The 
oscillation period is independent of (Iaf and solely determined 
by Al and the junction length L. The magnetic fields for 
dAF = 2.3nm and 5.8nm were divided by the geometry factor 
5 according to their shorter L. The side-minima are lifted up 
due to the measurement resolution of ~ 1 fiA. The calculated 
Ic{H) pattern (Eq. (3), dotted line) is offset by — 0.1/c//cO- 



tern for all interlayer thicknesses (Iaf matches the short 
J J model: 



Ic{H) = 1,0 



sm TT 



(3) 



with magnetic flux (f> = hqHL{2\l + d/ + dAp) (Lon- 
don penetration depth Al). This indicates (i) the flux 
penetrated cross-section area L(2Az, + + (Iaf) — 2LXl 
is unchanged for all magnetic fields, (ii) the Josephson 
coupling is uniform over the junction area, and (iii) the 
antifcrromagnctic Cr layer docs not modify the magnetic 
screening of the superconducting electrodes. Note that 
this is different for SI-F-S JJs, where flux focusing of the 
F-layer, the proximity- induced increase of Ai,, and reori- 
entation of the magnetic moments may lead to reduced 
oscillations periods for thicker dp. 



£,AF = 2.13 nm (set 2) and ^af = 3.44 nm (set 3). These 
coherence lengths are larger than in S-AF-S JJs with 
highly disordered FeMn interlayer, where S^af = 1.2 nm 
(after factor 2 definition correction of ^af) was observed 
[9]. The measured critical current density for JJs with- 
out chromium interlayer but the same tunnel barrier ox- 
idation conditions is roughly a factor 4 larger than the 
extrapolated jddAF = nm) from the sets with thinner 
dAF (set 1 and 2) , indicating additional scattering arising 
at the Cr interfaces. 



Specific resistance 

To check for a change in the chromium resistivity 
DC transport measurements were performed on planar 
chromium films on Si02 substrates. Fig. 3(b). The spe- 
cific resistance p of the AF-layer could not be directly 
determined from SI-AF-S junctions as the tunnel barrier 
resistance masks the serial resistance of the Cr layer. The 
residual resistivity ratio is nearly 1, indicating that p is 
set by the in-plane grain boundary scattering rather than 
by the temperature-dependent electron-phonon interac- 
tions. However, the current transport in SI-AF-S JJs 
is determined by the out-of-plane resistivity, for which 
grain boundary scattering is less important. The slope 
of jc{dAF) changed around 20 nm (set 2), probably due 
to modifications in p. The increase in ^af for thicker 
dAF, i.e. set 3, is consistent with the lower p(dAF) for 
these thicknesses. 

The estimation of the mean free path I in terms of 
the free-electron expression I = h/pe^ (v^Tr/n)^'''^ with 
carrier density n = 1.6 • 10^^ ni~^ yield £ = 2.2 nm, in- 
dicating that the dirty limit condition is not strictly ful- 
filled. Assuming a typical Fermi velocity for 3o? metals 
vp = 2-lO^m/s we obtain Z^AF — 14.6cm^/s. Taking the 
bulk E,^ yields ^af = 2 nm, which is fairly comparable 
to our ^AF determined from sets 2 and 3. This rough 
estimation docs not account for thin film modifications 
like interface-induced changes of the magnetic structure, 
interface diffusion, or domain formation, to name a few, 
which may vary Daf or Eex and hence ^af- 



B. Interlayer thickness dependence 

The jc{dAF) dependence could be measured over a 
range of four decades of jc, see Fig. 3(a). The IV- 
characteristics of all junctions were symmetric and repro- 
ducible after thermal cycling, thus flux trapping effects 
were unlikely to occur. The function 

jc ^ exp {-dAF/iAF) (4) 

[sec Eq. (2)] was fitted to the data. To avoid current dis- 
tribution effects only small JJs were considered, yielding 



C. Thin chromium interlayers 

JJs with thin Cr layers ranging from 2-6nm (set 1) have 
a spread of at least 30% in for the same dAF and a very 
weak or no thickness dependence, see inset of Fig. 3. The 
effective J J length of these junctions is ~ 2Aj, making 
long JJs effects unlikely. The spread of jc can be related 
to some variations of (i) the AF interlayer thickness, (ii) 
the magnetic structure, or (iii) the Josephson coupling. 
The normal state resistance i?„ was ^ 160 roD, for all 
these JJs, ruling out a variation of the junction area as 
a simple explanation. 
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FIG. 3: (Color online). jc{dAF) of 50^tm^ (set 1) and 500 ^tm^ 
(sets 2 and 3) JJs (a) and specific resistance p of single Cr films 
(b). The dotted line indicates the transition from the short to 
the intermediate JJs length regime. The inset depicts jcidAp) 
of set 1 on a linear scale. Dashed lines are fits to Eq. (4) for 
short JJs from sets 2 and 3. 



Variation of interlayer thickness 



The Cu layer at the bottom interface provides the 
growth of a uniform Cr layer. The strong spreads in jc 
require an interlayer thickness variation of roughly half 
the coherence length S,af, i-e. ~ Inm, and have not been 
observed in previous experiments with otherwise similar 
SI-F-S JJs (F=3d magnets) with a considerably shorter 
[G, 17, 18]. The steeper slope oi jcidAp) for thicker Cr 
layers makes jc even more prone to variations in (Iaf , but 
such large variations in were not seen for (Iaf > 6 nm. 
Furthermore, the magnetic diffraction pattern indicates 
uniform flux penetration in the barrier for all AF thick- 
nesses, see Fig. 2, i.e. uniform thickness of the Cr layer. 



Variation of the magnetic structure 



The appearance of some stochastically localized mag- 
netic states in the Cr (e.g. due to frustration) predomi- 
nantly for thin AF-layers would also affect the transport 
properties. The zero-field cool-down even facilitates the 
formation of such magnetic defects. A J J with a net 
magnetization in the interlayer deviates from the Ic{H) 
pattern Eq. (3), as observed in SI-F-S JJs [19]. As our 
junctions have a standard Ic{H) the local variations of 
the magnetic structure, if present at all, have to be small. 



Variation of Josephson coupling 

Considering the microscopic even/odd model the crit- 
ical current over a J J can be written as: 



Ic{Hx,Hy,T) 



jr(a:,y,T) sin Uo + ^ + 



H,x,x Hyy 



dx dy 



with junction length L, width w, phase 0o and j^'" being 
the critical current density of even or odd layers, respec- 
tively. For an odd number of AF layers the S-AF-S JJs 
ground state phase is a function of the temperature [4] , 
and at our measurement temperature of 4.2 K the odd 
layer parts could still be in the coupled ground state. 
Nevertheless, the jcS for even and odd number of AF 
layers should differ in magnitude. 

We assume that our SI-AF-S JJs have some atomic scale 
roughness in dAF- For similar absolute and an 
equal distribution of even and odd number of AF lay- 
ers and the odd AF layers mediating tt coupling, the 
Ic{Hx, Hy, T) curves would vanish at zero magnetic field 
due to local cancellation of the critical current. More 
generally, if the fraction k of the junction area has an 
odd number of AF layers mediating tt (0) coupling, the 
maximum critical current Ico is [(1 — =F K\j°\] Lw, 

the integral over the axis perpendicular to the applied 
magnetic field (assuming a uniform distribution of even 
and odd layers) yields the locally averaged critical cur- 
rent density (1 — n)\ic \ T >Ajc\; ^-nd a total phase ktt (if 
odd layers are tt coupled). Thus, the Ic{H) pattern still 
has the conventional form given by Eq. (3). 
Our integral supercurrent measurement determines the 
averaged critical current density and the Ic{H) curves of 
our samples agree well with the conventional diffraction 
pattern, see Fig. 2. Polycrystalline samples like ours 
with areas being much larger than the in-plane grain size 
(order of 10 nm in sputtered Cu/Cr samples [20]) seem 
not to be suitable for the verification of the even/odd 
layer model by measurements of Ic{H), Ic{T), or Ic{dAF)- 
The possible presence of an incommensurable SDW in Cr 
modifies the spin orientation and is an additional source 
for deviations from the even/odd model. 

Nevertheless, for the thinnest layers (set 1) the par- 
ity of the Cr layer number may set the local criti- 
cal current density due to atomic scale roughness and 
thereby stochastically reduce the maximum Jco, whereas 
for thicker layers (sets 2 and 3) the atomic layer model 
looses its validity. 



D. Antiferromagnetic order of chromium 

The antiferromagnetic order of the Cr-layer was indi- 
rectly checked by SQUID magnetometry measurements 
of unstructured reference SI-Cr(4-24 nm)/Fc(3.4 nm)-S 
films, being otherwise similar to the SI-AF-S Josephson 
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FIG. 4: (Color online). Magnetization of SI-Cr/Fe-S as a 
function of the in-plane magnetic field for several (Iaf at 15 K. 
Exchange bias He (stars) and coercivity field He (circles) are 
depicted in the inset. The paramagnetic electrodes tilt the Fe 
magnetization towards the field-axis. 



junctions, sec Fig. 4. In AF/F bilayers the unidirectional 
exchange bias effect induces for the F-layer (i) a shift of 
the magnetic hysteresis on the field axis by the exchange 
bias field Hf. and (ii) an enhancement of its coercive field 
He [11]. In most cases of non- itinerant and itinerant AFs 
the exchange bias field is in the opposite direction 
than the cooling field. A positive exchange bias has 
been observed on some samples, and is caused by the mi- 
crostructure and amplitude of cooling field [21]. In the 
case of Cr an oscillatory dependence of due to reori- 
entation of the SDW with thickness [ I _] and temperature 
[22] was reported. He depends strongly on the quality of 
samples and the properties of the interfaces, e.g. atomic- 
scale roughness inherent to sputtered films may weaken 
the exchange bias He- Easier experimentally detectable 
is Hcj which increases for thicker Cr thickness [23]. 

The samples were field cooled in 1 T from room tem- 
perature down to 15 K. Fig. 4 shows the widening of the 
hysteresis, i.e. i/c, with increasing Cr interlayer thick- 
ness and a small positive exchange bias He- The posi- 
tive sign of the exchange bias indicates antiferromagnetic 
coupling at the Cr/Fe interface. The saturation magneti- 
zation docs not increase for thicker Cr films and is solely 
contributed by the Fe film. We estimate an atomic mag- 
netic moment of ~ 2.2-2.3 /xb per atom, as expected for 
Fe atoms. The 4nm Cr layer still exerts an exchange bias 
on the Fe layer and the coercive field is He « 2.5 mT, be- 
ing comparable to single Fe films [24]. Polycrystallinc 
structure, spin leaking, and paramagnetic contributions 
of the thick electrodes soften the magnetically hard Fe 



layer. The boundary conditions of Cr modify the prop- 
agation direction and the spin orientation of the SDW, 
thus it may differ in SI-AF-S from the one in SI-AF/F- 
S samples. We note that the SI-AF/F-S samples were 
field cooled to verify the magnetic properties, whereas 
the zero-field cooled AF-JJs may show an enhanced AF- 
domain formation. 



IV. CONCLUSIONS 

In summary, transport studies of SI-AF-S junctions 
showed a uniform Josephson coupling, a very weak de- 
pendence of je on (Iap up to about 6 nm, and a decaying 
je for thicker (Iaf- The antiferromagnetic order of Cr 
was indirectly verified using SI-AF/F-S structures. Cr 
has a coherence length of S^Ap ~ 2.13-3.44 nm. 
A chromium dioxide, Cr02, interlayer between supercon- 
ducting electrodes was reported to mediate Josephson 
coupling over hundreds of nanometers [25]. The very 
weak damping of the critical current for thin Cr films 
{cLaf < 6 nm) indicates a weak singlet pair breaking ef- 
fect for nanometer thin clusters or layers of remaining 
metallic chromium in the Cr02. 

In future work on SI-AF-S samples the coupling of spin 
density (AF) and plasma (S) waves may be studied as 
the large subgap resistance facilitates dynamic transport 
studies. The magnetic field and temperature dependence 
of le should be measured in epitaxial S-AF-S samples and 
compared to different thickness regimes in order to verify 
the low-temperature anomalous behavior of le [3] . 
A thin layers of Cr may be used (i) to act as oxygen 
diffusion barrier while keeping a large superconducting 
gap close to the Cr interface, and (ii) for superconduct- 
ing spin valve structures, e.g. artificial antiferromagnets 
of F-AF-F type [24] sandwiched by S-clcctrodes as the 
damping of je in Cr is relatively weak. 
Critical current diffraction measurements of JJs are sen- 
sitive to magnetic remanence of the interlayer [1!J] and 
gave no indication of an intrinsic magnetic field due to 
the Cr layer. Both the relatively weak oxygen affinity 
and the exchange bias of chromium indicate that thin 
Cr films on top a superconducting metal may reduce the 
number of surface spins in SQUIDs [2()] and magneti- 
cally control the remaining spins at the interface at mK 
temperatures without changing the overall magnetic-field 
characteristics of the device. 
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